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ABSTRACT 

Long gamma-ray bursts (GRBs) have been suggested to occur preferentially in low metallicity environment. 
We discuss the possibility and theoretical aspects of using Lyman a emission properties of long GRB host 
galaxies as a metallicity indicator of high redshift GRB environments, where direct metallicity measurements 
are not easy. We propose to use the fraction of Lya emitters (LAEs) in long GRB host galaxies as a function 
of UV luminosity, which can be compared with star-formation-rate weighted LAE fraction of Lyman-break 
galaxies as the standard in the case of no metallicity dependence. There are two important effects of metallicity 
dependence of long GRB rate to change the LAE fraction of host galaxies. One is the enhancement of intrinsic 
Lya equivalent width (EW) by stronger ionizing UV luminosity of low metallicity stellar population, and the 
other is extinction by interstellar dust to change the observable EW. Based on a latest theoretical model of 
LAEs that reproduce observations, we argue that the latter is likely to work in the opposite direction to the 
former, i.e., to decrease LAE fraction if GRBs preferentially occur in low-metallicity environments, because of 
the clumpy inter-stellar medium effect. The high LAE fraction of GRB host galaxies indicated by observations 
is quantitatively be explained by the LAE model if GRBs occur when Z < 0.1Z Q , although this result is still 
indicative because of the limited statistics and theoretical uncertainties. This result demonstrates that the LAE 
statistics of GRB hosts may give us useful information in the future. 
Subject headings: gamma rays: bursts — galaxies: high-redshift 



1. INTRODUCTION 

Long gamma-ray bursts (GRBs; hereafter GRB means 
long gamma-ray burst unless particularly mentioned) are the 
brightest astronomical transient event, giving us an important 
laboratory of high energy astrophysics in extreme conditions, 
and an important tool to probe the high redshift universe. As- 
sociation of some observed GRBs with energetic type Ic su- 
pernovae (SNe) is considered to be an observational evidence 
that GRBs origi nate from core-collapses (CCs) of ver y mas- 
sive stars (e.g. iHjorfh et al. l2003HStanek et al. 120031) . How- 
ever, the occurrence rate of GRBs is much lower than that of 
normal CC SNe, and the condition required for a GRB to oc- 
cur from a CC still remains as one of the most outstanding 
questions about GRBs. 

Theoretical studies on possible G RB progenitors and pro- 
ductio n mechanisms of GR B s (e.g. MacFadven & Wooslev 
1999t lYoon&Langerl 12005b lYoon et al. I 120061: 
Woosley & Hergerl 120061) suggest that low metallicity 
may be a necessary condition for the progenitor to produce 
a GRB event. It has also been suggested from observations 
that GRB host galaxies are systematically fainter than those 
in the case that GRBs are unbiased tracer of star formation or 
CC SN rate, indicating that GRBs may preferentially occur 
in low metallicity environment because fainter or smaller 
galaxies generally have lower metallicity (Le Floch et al. 
2003; Fruchter et al. 2006; Wolf & Podsiadlowski 2007). 
These interpretations have also been supported by studies 
using theoretical models of galaxy formation and evolution 
(Nuza et al. 2007; Lapi et al. 2008). Furthermore, Stanek et 
al. (2006) reported that the distribution of metallicity of five 
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GRB host galaxies at z < 0.25 is significantly biased towards 
low metallicity compared with the expectation when GRBs 
are unbiased star formation tracer. 

However, spectroscopic estimates of metallicity are avail- 
able only for galaxies at low redshifts (z < 1, Stanek et al. 
2006; Savaglio, Glazebrook & Le Borgne 2009), while ma- 
jority of GRBs occur at higher redshift. At higher redshifts, 
the metallicity of host galaxies is discussed based on host 
galaxy luminosity or mass, with the empirical relation of 
luminosity-metallicity and mass-metall icity (L-Z and M-Z, 
iFruchter et alTI l2006t IChen et al. l |2Q09l). However, there is 
a large scatter in the L-Z relation. Although t he M-Z relation 
of ga l axies is tighter t han the L-Z relation (Tremon ti et al. I 
2004; Erb et al. 2006), rest-frame infrared luminosity must 
be measured to reliably estimate galaxy stellar mass, which 
requires mid-infrared observations such as those by Spitzer. 
Therefore it would be useful if there is another and indepen- 
dent metallicity-sensitive property of galaxies that can rela- 
tively easily be observed even at high redshifts (z > 1) where 
the majority of GRBs occur. In this paper we consider to use 
the Lya emission of galaxies as such an indicator. 

Lyman alpha emitters (LAEs) can be detected from 
very large redshif t s by their strong Lya emission lines 
( Ouch i et al. 1120081: iMartin et al. 112008b ISawicki et ani2008t 
i Finkelstein et al. 1 12009b iNilsson et al. 1 120091: IShioya et all 
2009, and reference therein), including the currently high- 
est redshift galaxy confirmed spectroscopically at z = 6.96 
(lye et al. 2006, Ota et al. 2008). LAEs [typically defined 
as galaxies having rest-frame equivalent width (EW) > 20 
A] ar e found only in metal poor galaxies in the local uni- 
verse (Chariot & Fall 111993b . Studies of spectral energy dis- 
tribution (SED) of LAEs at z > 3 have suggested that they 
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are typically younger and less-massive than UV-continuum 
selected galaxies [e.g., Lyman-break galaxies (LBGs)1 at 
simila r redshifts (e.g. [pawiser et al. 200 61: iNilsson et al71 
I2007t iFinkelstein et al. I l2007t lLai et al. 1 120081) . indicating 
low metallicity of LAEs. Since low metallicity stars emit 
more ionizing UV photons, it is theoretically expected that 
low metallicity galaxies tend to have large Lya EW. Metallic- 
ity may affect the Ly a pro perties also through extinction by 
interstellar dust (see 32.31 ). Therefore, although some other 
effects such as dynamics of inter-stellar medium (ISM) may 
well play an imp ortant role to determine the Lya properties 
(lAtek e t al. 2008, and references therein), the probability of 
GRB host galaxies being LAEs could be an indicator of metal- 
licity dependence of GRB rate that can be used even at high 
redshifts. In fact, observations have indicated that the LAE 
fraction in GRB host galaxies may significantly be higher than 
that of field galaxies. Though the sample is still small, almost 
all of high redshift (z > 2) GRB hos ts are LAEs dFvnbo et al. I 
|2002j,[2603; Jak obsson et aL~l l2005). while only about 10% - 
25% of L BGs at similar redshifts have such strong Lya emis- 
sion lines (Sha plev et al. 1120031: iReddv et al. 112008b . 

We will discuss the theoretical aspects of Lya emission of 
GRB host galaxies, and prospect for using it as a tool to get 
information about GRB progenitor from future observations, 
especially about the metallicity dependence, based on the lat- 
est developments of observations and theoretical models of 
high redshift LAEs. Lapi et al. (2008) also discussed LAEs 
as GRB host galaxies based on the LAE model of Mao et al. 
(2007), but their model assumes that 'intrinsic' ionizing pho- 
ton luminosity (i.e., that originally from stars without taking 
into account extinction or absorption) in each galaxy is sim- 
ply proportional to star formation rate (SFR) in the galaxy, 
without dependence on age or metallicity. It is difficult to 
examine the dependence of LAE fraction on metallicity with 
such a modeling. Here we use one of the latest models of LAE 
luminosity functions (LFs) of Kobayashi, Totani, Nagashima 
(2007; 2009, hereafter KTN07 and KTN09). This is based 
on a hierarchical clustering model of galaxy formation in the 
framework of cosmological structure formation theory, and 
the intrinsic ionizing luminosity is calculated using star for- 
mation history taking into account metallicity evolution (see 
32.21 ), and hence dependence of Lya luminosity on properties 
of stellar population can be examined in a realistic way. 

In $2] we discuss general theoretical aspects of Lya emis- 
sion and its dependence on galaxy properties and introduce 
the LAE model used in this work. In |3j we discuss obser- 
vational strategies to extract useful information about GRB 
progenitors from future data sets, with least dependence on 
theoretical uncertainties. In we quantitatively investigate 
how the fraction of LAEs in GRB host galaxies could be af- 
fected if the GRB rate is dependent on metallicity. In <(5] we 
compare our model prediction to current observed data set, 
and also discuss the age effect of the GRB progenitor and 
the selection effect on host galaxies induced by extinction of 
GRB optical afterglows. We will summarize our conclusion in 
§6\ We assume ACDM cosmology with £1 A = 0.7, S1 M = 0.3, 
and Ho = 70 km s" 1 Mpc" 1 , and use the AB magnitude system 
throughout this paper. 

2. THEORETICAL BACKGROUND AND THE MODEL OF LAEs 

2.1. The Galaxy Formation Model 

To construct a mock numerical catalog of galaxies, we 
utilize one of the latest hierarchical clustering models 



(so-called semi-analytic m odels) of galaxy formation by 
Nagashima & Yoshii (2004). This model, which we call the 
Mitaka model here, produces numerical catalogs of galaxies 
and their evolution by semi-analytic computation of merger 
history of dark matter (DM) halos and phenomenological 
model calculation of evolution of baryons. The computation 
of DM halo merger history is based on the standard theory 
of structure formation driven by cold dark matter. The model 
calculation of baryon evolution within DM halos includes ra- 
diative cooling, star formation, SN feedback and galaxy merg- 
ers. Metal enrichment history (chemical evolution) of inter- 
stellar medium is calculated self-consistently by the star for- 
mation history in each galaxy. 

The Mitaka model can rep roduce wide variety of obser ved 
properties of local galaxies (Nagashima & Yoshi il 120041) a s 
well as those of high redshift LBGs ( Kashika wa et al. II2006I) . 
In Fig. Q] we show cosmic star formation history (CSFH) 
computed in the Mitaka model with and without metallicity 
limit. The metallicity limited CSFHs (Z < 0.4 and O.2Z ) 
have their peak at higher redshifts than that of total CSFH. 
The evolution of cosmic SFR and dependence on the metal- 
licity limit are similar to those in othe r models utilized in 
studies of GRB rate evolution (e.g. iDaigne et al. 1 120061 : 
lLanger& Norman 1 120061: lLapi et al. 1 120081) indicating that 
our model CSFH is broadly consistent with the GRB flux or 
redshift distributions, within the uncertainties about GRB lu- 
minosity function. 

2.2. Lya Photon Production 

KTN07 extended this model to describe LAEs. Most of 
previous studies simply assumed that Lya luminosity is sim- 
ply proportional to SFR, but the KTN07 model incorporates 
the effect of dust extinction with amount different from that 
for UV continuum, and the effect of outflow. The KTN07 
succeeded to reproduce the Lya LF at various redshifts. Fur- 
thermore, KTN09 compared this model comprehensively with 
various observed data of UV continuum LF and equivalent 
width distribution of LAEs, and found a good agreement. We 
use the latest KTN09 model for the analysis in this paper, and 
see this paper for detailed description of the LAE modeling. 

There are several important ingredients to determine the 
Lya luminosity from a galaxy. The first is the ionization lumi- 
nosity (A < 912 A) that is calculated by star formation history 
and metallicity in each galaxy in the KTN09 model. Then the 
model predicts the intrinsic Lya luminosity assuming that a 
part of ionizing photons are absorbed within the galaxy and 
reprocessed into Lya photons by the case B recombination. 
The fraction of ionizing photons that are absorbed by neutral 
hydrogen to produce Lya is treated as a p art of the overall 
normalization factor of Lya (/o, see 32.31 ). The 'intrinsic' 
EW is then determined by the ratio of the intrinsic Lya lu- 
minosity to the UV continuum luminosity. Fig. |2] shows the 
evolution of the intrinsic EW for single burst stellar popula- 
tion or the case of constant star formati on rate fo r a var iety 
of metallicity, calculated by the model of Schaerer] d2003l) as- 
suming the Salpeter IMF in the range of 1-100 M@. Since 
stellar populations with low metallicity emit more ionizing 
photons compared with UV continuum photons around Lya, 
EW becomes larger. This effect is significant in young stellar 
populations (< several xlO 6 yr). LAEs having such young 
age e s timates have ac ually been observed (e.g. ICharv et al. I 
120051: lLai et al. 1 120071) . and they are important to reproduce 
large EW LAEs in the model of KTN09. This is one of the 
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important motivations to consider LAEs as the metallicity in- 
dicator of GRB progenitors. 

A caveat here is that the stellar evolution theory predict- 
ing the ionization luminosity may still be highly uncertain at 
such young ages, because of e.g., the treatment of the atmo- 
sphere of Wolf-Rayet sta rs (Garcia- Vargas, B ressan & Diaz 
1995) and stellar rotation ( Vaz que et al. 112007b . Our model is 
based on the widely used model of Schaerer (2003), and qual- 
itative trend that younger and lower metallicity stellar popula- 
tion should have large intrinsic Lya EW is generally accepted. 

2.3. Lya Escape Fraction and Extinction Effect 

The observable Lya luminosity is determined by the intrin- 
sic luminosity and the escape fraction of Lya photons from 
the galaxy, and KTN09 considered the effects of extinction 
by dust. KTN09 found that the extinction effect is especially 
important to successfully reproduce various observations of 
LAEs. In KTN09 model, the Lya luminosity of a galaxy in- 
cluding the extinction effect is computed as: 

- l-exp(-r Lya ) rmax 

^Lya-J0 L Lya (I) 

where TL yQ is an optical depth of the galaxy for Lya photons, 
and L™* is Lya emitted in the galaxy if all ionizing photons 
are reprocessed into Lya photons. The parameter T\^ a is pro- 
portional to that of UV continuum, Tuv, and both scales with 
the metal column density of galaxies calculated in the Mitaka 
model. The overall normalization factor (fo) and normaliza- 
tion of r Lya are determined by fitting to the observed LAE 
Lya luminosity functions. 

Because of the resonant scattering by neutral hydrogens, 
Lya photons would take vastly different paths in a galaxy 
from those of UV continuum photons, and hence TL yct can 
be very different from tuv- It is theoretically highly uncertain 
whether this effect reduces or enhances Lya EW, and there 
are two extreme possibilities: extinction of Lya photons can 
be much more significant (in the case of homogeneous ISM, 
Ferland & Netzer 1979) or much less (in the case of clumpy 
ISM, Neufeld 1991; Hansen & Oh 2006) than that of UV 
continuum. Recent observations by Finkelstein et al. (2007, 
2008, 2009) have indicated that the clumpy ISM effect is in 
fact working in high-z LAEs. This has theoretically been sup- 
ported by KTN09 who independently found that the clumpy 
ISM case (i.e. TL yQ < tuv) is favored to reproduce the exis- 
tence of large EW (> 240 A) LAEs, EW-reddening correla- 
tion, and EW-Luv correlation. If this interpretation is correct, 
we expect that galaxies with large metallicity (and hence a 
large amount of dust) tend to have larger Lya EW, which is 
in the inverse direction to the effect of ionization luminosity 
mentioned above. 

It has been argued that LAEs are less dusty galaxies, both 
theoretical ly (e.g. iFerl and & Ne tzer 1 11979 ) and ob serva- 
tionallv dGawiser et al. 1 120061 120071: iGronwall et al 1 120071: 
iNilsson et al. I 12007b lLai et al. I 120081: lOuchi et al. 1 120081) . 
These claims appear in contradiction with the preference of 
the clumpy ISM dust suggested by recent studies of Finkel- 
stein et al. (2007, 2008, 2009) and KTN09. However, this can 
be understood by Fig. 7 of KTN09 where the plot of extinc- 
tion (A1500 at 1500 A) versus EW is shown for the model and 
the observed data. If one selects LAEs with modest values of 
EW (< 100 A), the extinction is relatively small, while the 
observed data of Finkelstein et al. (2009) indicate that LAEs 



with large EW (> 100 A) are more dusty, which is also re- 
produced by the KTN09 model. Since LAEs with EW < 100 
A are more abundant, a conclusion that LAEs are not dusty 
could be derived depending on the sample size and selection 
criteria of LAEs. 

2.4. Remarks about LAE Theory in the Context of This Work 

As discussed above, the consequence of the metallicity de- 
pendence of GRBs on the Lya properties of host galaxies 
would be rather complicated and theoretically still uncertain. 
We use the KTN09 as the guideline in this work, but it should 
be noted that there are still large uncertainties in understand- 
ing of LAEs. Some of the theoretical predictions may be for 
particular cases, or be changed by future studies. However, 
our result would still be useful as the first step to consider fu- 
ture possible use of LAEs to derive some information about 
GRBs. We also emphasize that the KTN09 model is one of 
the latest ones, and unique to reproduce observed data consis- 
tently at various redshifts for all of the LAE LFs in Lya and 
UV continuum luminosities, EW distributions, and Luv ver- 
sus EW correlations, taking into account the selection criteria 
of LAEs in each observation. 

3. OBSERVATIONAL STRATEGY 

It is not straightforward to extract useful information from 
the statistics of LAE fraction in GRB host galaxies. Even if 
GRBs are simply tracing SFR, the LAE fraction for a given 
sample of GRB host galaxies is weighted by SFR of galaxies, 
and hence different from that for a sample of galaxies found 
by flux-limited galaxy surveys. If we have another transient 
source that traces SFR, such as CC SNe, we can compare host 
galaxies of GRBs and the another source. However, LAEs are 
found at z > 2 and CC SNe cannot be detected at such high 
redshifts. 

We propose to compare the LAE fraction of GRB host 
galaxies with the SFR- weighted LAE fraction of high-redshift 
galaxies found by flux-limited galaxy surveys. It should be 
possible to construct a flux-limited sample of GRB host galax- 
ies, with the same flux limit as that of a galaxy survey sample. 
Then we expect that the LAE fraction of GRB host galaxies 
should be the same as the SFR-weighted LAE fraction of the 
galaxy survey sample, provided that GRBs trace SFR. In other 
words, a significant difference between the two indicates that 
GRBs do not simply follow SFR in galaxies. A caveat is pos- 
sible selection effects in detecting GRB host galaxies. If the 
detection efficiency of GRB hosts does not depend on galac- 
tic properties except for the flux limit for galaxies, the com- 
parison will not be biased. However, even if the prompt and 
afterglow emission properties do not depend on host proper- 
ties, the detection probability of optical afterglow could be 
affected by extincti on in galaxies. We will discuss this issue 
later in this paper ( {35.31 l. 

We need estimates of SFR of galaxies in the galaxy survey 
sample. One possibility is to use rest-frame UV luminosity as 
the SFR indicator, because this is always available for galax- 
ies found by the Lyman-break method (i.e., LBGs). However, 
UV luminosity may be affected by extinction, and a more re- 
liable SFR estimate is favorable. Such an estimate is also pos- 
sible, e.g., by SED fittings to multi-band photometric d ata or 
spectroscopic observations (e.g. IShapley et al. 1 120051) . The 
number of galaxy data sets with such rich information is in- 
creasing thanks to the recent developments of high-redshift 
galaxy searches. It should be noted that an observational SFR 
estimate is generally an average of SFR until the time of the 
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observation on a time scale that is dependent on SFR indica- 
tors. Here, we implicitly assumed that the evolutionary time 
scale of SFR of a galaxy is longer than the time scale of SFR 
indicator, so that the SFR estimate is not bias ed by star for- 
mation history. This issue will be discussed in jj5.2| 

The upper panel of Fig. [3] shows the simple LAE fraction 
in number and SFR-weighted LAE fraction of galaxies in the 
Mitaka model at z= 3, as a function of restframe UV (A= 1500 
A) magnitude, Muv- Two different threshold values of Lya 
EW are adopted as the definition of LAEs. It can be seen that 
the SFR-weighted LAE fraction is significantly higher than 
the simple number fraction at a fixed value of Muv- This 
can be understood by the effect of clumpy dust for the ex- 
tinction of Lya photons. For a fixed Muv, galaxies with large 
extinction should have larger (extinction-corrected) SFR, and 
large extinction increases Lya EW by the clumpy dust effect. 
Therefore, even if there is no metallicity effect about GRB 
rates, and GRBs simply trace SFR, we still expect a higher 
LAE fraction of GRB host galaxies than the simple number 
fraction of LAEs in galaxy surveys. This result also implies 
that we should not use simply Muv as SFR indicator; if Muv 
is a good SFR indicator, we expect the same value for the 
SFR-weighted and simple number fractions at a fixed Muv- 

The lower panel of Fig. [3] shows the UV continuum LF of 
galaxies in the Mitaka model. The UV-luminosity-weighted 
LF and SFR-weighted LF are also shown. The luminosity dis- 
tribution of GRB host galaxies should be the same as the SFR- 
weighted LF if GRBs simply trace SFR. The SFR-weighted 
LF shows some deficit in the largest UV luminosity range, 
compared with the UV-weighted LF. This can be understood 
by the fact that galaxies with large extinction tend to have 
smaller observed (i.e., extinction-uncorrected) UV luminos- 
ity, and hence higher SFR compared with UV luminosity. A 
comprehensive comparison between the model prediction of 
UV continuum LF in the Mitaka model with observations are 
presented in KTN09. 

4. LOW METALLICITY PREFERENCE OF GRBs AND LAE 
FRACTION 

Now we consider the effect of the metallicity dependence 
of GRBs on LAE fraction. As the first simple model, We 
assume that GRBs occur with a rate proportional to SFR but 
only when the ISM metallicity of the galaxy is smaller than 
a critical value, and GRBs do not occur otherwise (the step 
function model): 

* GRBOC {o FR ' Z>Ze!-|- (2) 

This simple model of GRB rate was tested in number of stud- 
ies and reproduced various statistical properties of GRBs such 
as redshift distribution and logV - logP plot within the range 
of uncertainty when it is combined with realistic model of star 
formation history and metallicity of galaxies (e.g. lLapi et al~l 
l2008h . 

Our interest is in the fraction of LAEs in GRB host galax- 
ies, and hence we do not need to consider the normalization 
(the proportionality constant) of Rgrb, luminosity function or 
spectrum of GRB prompt emission. Here we calculate SFR 
by that of the model galaxy at the time of an observation, and 
implicitly assume that the SFR evolution time scale is longer 
than the typical age of the GRB progenitor. If SFR is chang- 
ing with a shorter time scale than the progenitor age, the real 
GRB rate could be different from that cal cula ted by this for- 
mulation. This issue will be discussed in 35.21 



We then calculate the fraction of LAEs in GRB host galax- 
ies with various values of Z cl -i t using the Mitaka model at z = 3, 
and the result is shown in Fig. |4]as a function of UV luminos- 
ity, for the two different threshold values of LAE EW as in 
Fig. [3] The LAE fraction without the metallicity dependence 
(Z CI it = oo) is plotted together. Compared with the case of no 
metallicity dependence, the predicted LAE fraction becomes 
higher at large UV luminosity (Muv 19) and lower at faint 
(Muv ^ -19). These trends can be understood as follows. 

Galaxies with relatively small amount of extinction are 
dominant in the large UV luminosity range, since large extinc- 
tion would significantly reduce the UV luminosity. In such a 
regime, the effect of extinction by clumpy ISM on Lya EWs is 
not significant, and EWs are mainly determined by the stellar 
population, i.e., larger EWs for younger and lower metallic- 
ity galaxies. Therefore, if GRBs trace low metallicity stellar 
populations, an enhancement of Lya EW is expected. On the 
other hand, significantly extincted galaxies populate the low 
UV luminosity range, and the effect of clumpy ISM dust on 
Lya EW is important. The dust/gas ratio should be larger 
for large metallicity galaxies , and hence our model predicts 
large Lya EW for metal-rich galaxies in this regime. 

Therefore, the key is that UV-faint galaxies have larger ex- 
tinction on average, combined with the clumpy ISM effect. 
This may sound strange for the reader, because it is known 
that galaxies with larger stellar mass tend to have larger ex- 
tinction. However, we emphasize that the UV luminosity dis- 
cussed here is the apparent luminosity that is not corrected 
for dust extinction. Galaxies with larger stellar mass does not 
necessarily have larger UV luminosity, due to the larger dust 
extinction in UV bands. To summarize, the two competing ef- 
fects of metallicity on the LAE fraction discussed in §|2] i.e., 
stellar population and extinction, appear in different ranges of 
UV luminosity of host galaxies. 

We also calculate the LAE fraction of GRB host galaxies 
integrated by UV magnitude in the range corresponding to 
the observed data with the threshold EW of 10 A, and show 
in the left panel of Fig. [5] as a function of the critical metal- 
licity. The LAE fraction rapidly increase with decreasing Z cl i t 
atZ cr i t ~ 0.1Z Q . This metallicity scale is interestingly similar 
to those suggested by theoretical studies of stella r evolution 
dWooslev & Hergerll2006t lYoon & Langer 1 120051) and spec- 
troscopically measured metallicity of loca l GRB host galaxies 
dStanek et al. ll2006t[Savaglio et al. 112009b . This sensitivity of 
the LAE fraction for the metallicity scale of our concern in- 
dicates an interesting opportunity to constrain the GRB pro- 
genitor by future observations of high redshift galaxies and 
GRBs. 

It should be noted that the behavior of LAE fraction as 
a function of Z alt is not simple, but it depends on the EW 
threshold and UV magnitude range of the galaxy sample con- 
structed. As discussed previously in this section, LAE frac- 
tion increases with decreasing Z cr ; t in the bright UV magni- 
tude range, and decreases in the faint UV range. Hence, if we 
consider only brightest galaxies, the LAE fraction converges 
to unity due to the stellar population effect, which is the case 
of the left panel of Fig. [5] (The effect of clumpy ISM which is 
into the opposite direction to the ionization luminosity effect 
is also seen in the range of log Z a i t ~ -0.6 - -0.3.) However, 
if we take fainter UV magnitude range into account, the be- 
havior of LAE fraction becomes complicated. LAE fraction 

1 The amount of interstellar dust is assumed to be proportional to the gas 
metallicity in the Mitaka model. 
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increases with decreasing Z cl ; t at large Z cl ; t , but the trend turns 
over and it decreases with decreasing Z cr i t at small Zdt, for 
some values of the EW threshold, as shown in the right panel 
of Fig. [5] This point should be kept in mind when the fu- 
ture observations are used to constrain the nature of the GRB 
progenitor. 

The step function with the critical value of Z cr j t is likely to 
be too simple to describe the realistic GRB production effi- 
ciency as a function of metallicity. Therefore we also try a 
more smoothly varing function (the Gaussian model): 



Rcrb exp ( -^2^ x SFR, 



(3) 



and the LAE fraction with this model is plotted against <7z 
instead of Z a - lt in Fig. [5] Qualitative result is not changed 
from the case of the step function model. 

5. DISCUSSION 
5.1. Comparison to Observations 

Here we compare our result with observations. Currently 
only one published result is available for the LAE statistics of 
GRB host galaxies (Jakobsson et al. 2005), which is based 
on seven GRB hosts with Lya observations including those 
of Fynbo et al. (2002, 2003). Their result is summarized in 
Table Q] In these seven host galaxies, five have strong Lya 
emission (EW > 10 A). Though Lya emissions were not de- 
tected in two host galaxies (000301c and 020124), the upper 
limits are not severe and we cannot exclude the possibility 
that these two are also LAEs having EW > 10 A. It is difficult 
to set strong constraints on Lya EW for galaxies fainter than 
^28 mag in optical (rest UV) bands by the current observing 
facilities, and hence we do not include the three faintest host 
galaxies in Table Q] (000301c, 020124 and 030323) in the fol- 
lowing discussion. Now we have four GRB host galaxies at 
Z ~ 2-3 with known Lya EWs of 14 A - 71 A. Absolute rest 
UV magnitudes of these are Muv ~ _ 22 - -21. It is obvious 
that the statistics is not yet sufficient, but here we compare 
these results with our model at z = 3 as the first trial. 

The observational constraints are indicated by shaded re- 
gions in Fig. |U using the confidence limits for small number 
statistics given by Gehrels (1986). In these plots adopting EW 
thresholds of 40 and 120 A, our model predictions are roughly 
consistent with the observed data, giving no constraint on Z cr j t . 
However, if we take the EW threshold of 10 A and the thresh- 
old absolute rest UV magnitude of Muv - 5 logh < -21 to 
match the sample of Jakobsson et al. (2005), Z a \ t < 0.1Z Q 
seems to be favored in our model (Fig. 0, which is interest- 
ingly consistent with other theoretical or observational indica- 
tions. However, the statistics is obviously not sufficient, and 
we need to await future observations to derive more reliable 
conclusions. 

5.2. Progenitor Age Effect on LAE Fraction 

As mentioned before, we have implicitly assumed that the 
evolutionary time scale of SFR in galaxies is longer than the 
age of the GRB progenitor and the age of stellar popula- 
tion used to observationally estimate SFR. However, if SFR 
is rapidly changing with a comparable or shorter time scales 
than these, our argument above could be affected. Since we 
expect strong Lya emission from very young stellar popula- 
tion having large ionizing luminosity (Fig. |2), the evolution- 
ary time scale of such galaxies might become as short as the 



time scales of the GRB progenitor or SFR indicators. 

Fig. 6 of KTN09 shows the distribution of mean stellar ages 
of LAEs in the Mitaka model, where the age is calculated by 
star formation history of each galaxy with a weight by ioniza- 
tion luminosity. Since ionization luminosity comes from mas- 
sive stars having short lifetimes, the ionization-luminosity- 
weighted age is about 10 6 6 yr for stellar population having 
a constant SFR with a duration much longer this time scale. 
According to Fig. 6 of KTN09, LAEs with large EW have 
significantly shorter age than this, indicating that their SFR is 
not constant on this time scale. 

If the GRB progenitors have a longer time scale to evolve 
into a GRB than this scale, the age effect is not negligible. 
According to the main-sequence lifetimes of massive stars 
(Schaller et al. 1992; Lejeune & Schaerer 2001), it may be 
the case. In galaxies whose mean stellar age is significantly 
shorter than the age of GRB progenitor, we expect a smaller 
GRB rate than simply expected from SFR because the progen- 
itors do not have enough time to produce GRBs, while such 
young galaxies have a high probability of being LAEs. This 
effect would then lead to apparent decrease of LAE fraction 
of GRB host galaxies, and hence it might be confused with 
the effect of the low metallicity preference of GRBs com- 
bined with the extinction by dust in the clumpy ISM. This 
is an interesting possibility, and should be kept in mind when 
we analyze future data sets of GRB host galaxies, though a 
quantitative prediction is beyond the scope of this work. 

The typical lifetime of stars contributing to the continuum 
UV luminosity at 1500-2800 A, which is a popularly used 
SFR indicator, is ~ 20 Myr (corresponding to ~ 10M Q stars, 
Madau, Pozzetti & Dickinson 1998). Therefore, it is also pos- 
sible that SFR of LAEs is changing with a time scale shorter 
than the time scale of the SFR indicator. One must be care- 
ful about this when the SFR-weighted LAE fraction of field 
galaxies is calculated as the reference value corresponding to 
the case of the SFR-tracing GRB rate. 

5.3. Dust Extinction of Optical Afterglows 

In most cases, detection of an optical afterglow is required 
to identify the host galaxy of a GRB at high redshift, and 
hence a sample of identified GRB host galaxies may be biased 
against dusty host galaxies in which afterglow light is signif- 
icantly attenuated. Since the extinction should have a signifi- 
cant effect on the Lya EW of host galaxies as discussed in $2] 
the extinction bias of GRB host galaxies could affect the LAE 
statistics of them. 

To investigate this effect, we employ a simple method simi- 
lar to that used to investigate the low metallicity preference of 
GRB events. We assume that we can identify the host galaxy 
of a GRB only when the extinction of the galaxy is less than 
a critical value, and hence the effective GRB rate becomes: 



/SFR, A v <A Vc 
I ' 



A v >Ai 



(4) 



The LAE fraction of GRB host galaxies at z = 3 predicted 
by this model is shown in Fig. [6] and the LAE fraction 
decreases as the afterglow extinction effect becomes strong 
(smaller Av.ait)- This is simply a manifestation of the effect 
of extinction by clumpy ISM dust on LAEs, i.e., enhancement 
of LAE EW by extinction, which is assumed in the model of 
KTN09. The change of the LAE fraction is into the same di- 
rection in all the magnitude range, and it can in principle be 
distinguishable from the effect of low metallicity preference 
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of GRBs showing increase and decrease of the LAE fraction 
with decreasing Z cr i t at Muv - 51og/z < — 19 and > — 19, re- 
s pectively (Fig. p}. 

Fyn bo et al. 1 12003) pointed out that the bias against dusty 
galaxies by the extinction of optical afterglows may be the 
reason of the observed high LAE fraction of GRB host galax- 
ies, because it has often been argued that LAEs are less dusty 
galaxies. However, as discussed in $2] recent studies have 
suggested that the effect of extinction on LAE EW is rather 
inverse and to increase EW by extinction in clumpy ISM. 
Therefore, if the indication of the clumpy ISM extinction is 
correct, the extinction of GRB afterglows is unlikely to cause 
large LAE fraction among host galaxies. 

6. CONCLUSIONS 

We discussed the theoretical aspects of using Lya emission 
properties of GRB host galaxies to get implications for the 
GRB progenitor. The fraction of LAEs in a sample of GRB 
host galaxies can be compared with the SFR-weighted LAE 
fraction of a sample of LBGs at similar redshift with a simi- 
lar magnitude limit, and a significant difference, if observed, 
would indicate that GRBs do not simply follow SFR and a 
new physical parameter, such as metallicity, is hidden in the 
relation between GRB rate and SFR. 

To make quantitative predictions we used one of the latest 
theoretical model of LAEs in the framework of hierarchical 
galaxy formation, which is successful to reproduce a variety 
of recent observed data about LAE statistics. We then tested 
the case that GRBs have preference to low metallicity envi- 
ronments. We found that there are two effects of the metallic- 
ity dependence on the LAE fraction of host galaxies. One is 
the intrinsic EW enhancement by low metallicity; low metal- 
licity stellar populations have larger ionizing UV luminosity 
and hence larger Lya EWs. The other is the extinction by 
ISM dust; low metallicity galaxies are expected to be less 
dusty and dust affects LAE EWs. According to the KTN09 
model as well as recent observations, Lya EWs may be en- 
hanced because of extinction by dust in clumpy ISM. There- 
fore, the low metallicity preference of GRBs would decrease 
the LAE fraction by this effect, which is into the opposite 
direction to that of the intrinsic EW enhancement. We pre- 
dicts that the intrinsic EW enhancement is more significant 
for galaxies with large UV luminosities, while the extinction 



effect dominates in UV-faint galaxies, separated by the border 
of Mu v - 5 log /; ~ - 1 9 at z = 3 . 

The observational indication of higher LAE fraction for 
GRB host galaxies reported by Fynbo et al. (2002, 2003) and 
Jakobsson et al. (2005) can be explained by our model by the 
intrinsic EW enhancement effect, if GRBs occur only in low 
metallicity environment of Z < Z cl ; t = O.1Z . However, the 
statistics is still poor and further systematic observations of 
GRB host galaxies are highly desired. 

We have also discussed the selection effect of GRB host 
galaxies by extinction of GRB afterglow flux, which is an- 
other possible explanation of large LAE fraction. However, 
recent observations and theoretical modeling indicate that the 
effect of extinction is to enhance Lya EW by clumpy ISM, 
and in this case the selection effect will decrease the LAE 
fraction. 

It should be noted that the theory of LAEs is still highly un- 
certain, and the quantitative predictions by the KTN09 model 
suffers from such uncertainties. However, the quantitative be- 
haviors about LAE emission of GRB hosts presented here 
based on one of the latest models of LAEs would be useful 
to start consideration about using LAEs to derive information 
about GRBs by the future data set. Observations of high red- 
shift LAEs are rapidly developing, and near-future observa- 
tions of LAEs will provide us with further knowledge about 
this galaxy population. We tried to outline the important theo- 
retical effects and aspects that should be considered when one 
uses Lya properties of GRB host galaxies to get any implica- 
tions about GRBs from future observations. 
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TABLE 1 

LONG GRB HOST GALAXIES WITH LYO EW MEASUREMENTS 



GRB 



redshift host mag Fi^ ya (erg s 1 cm 



) L Lya (evgs-') Lya EW rcsl (A) 



971214 


3.42 


R 


= 26.5 


~ 4.5 x 10" 


17 


~ 5.0 x 10 42 


~ 14 


000301c 


2.04 


K 


= 28.0 


< 3.5 x 10" 


17 


< 1.1 x 10 42 


< 150 


000926 


2.04 


V 


= 24.9 


1.6 x 10- 16 




5.1 x 10 42 


71 


011211 


2.14 


R 


= 24.9 


2.8 x 10~ 17 




1.0 x 10 42 


21 


020124 


3.20 


V 


= 29.3 


< 1.5 x 10" 


17 


< 1.4 x 10 42 


<22 


021004 


2.34 


B 


= 24.4 


2.5 x 10- 16 




1.1 x 10 43 


68 


030323 


2.14 


V 


= 28.0 


1.2 x ur 17 




~ 1.3 x 10 42 


~ 145 



NOTE. — Host magnitude of GRB 020124 is from Bloom et al. (2002), and other values are from Jakobsson et al. 12005, see their Table 4 and references 
therein). The magnitudes have been converted to AB magnitudes. 



DC 
LL 



10" 



total 
Z < 0.4 Z ffl 




3 4 
redshift 



FIG. 1. — CSFH computed in the Mitaka model. CSFH of all galaxies, and those of low-metallicity galaxies (Z < 0.4Zq and 0.2Zq) are shown by solid, 
dashed, and dotted lines, respectively. 
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Z = (Pop III) 
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log 10 [t/yr] 

FIG. 2. — The time evolution of intrinsic Lyce equivalent width (EW[^ a ) of stellar populations calculated by the stellar population synthesis model o flSchaererl 

12003) with the Salpeter IMF in the range of 1—100 Mq. The solid and dashed curves represent the EWj? a evolution for instantaneous starburst and constant 

star formation, respectively. Several curves are shown corresponding to different stellar metallicities: Z = 0.(1/2000)Zq , (1/50)Zq , (1/20)Zq , (Y/5)Zq , and 
Zq, from top to bottom. In this figure, all the ionization photons are assumed to be absorbed by neutral hydrogens to produce hya photons by the case B 
recombination. 




FIG. 3 . — Upper Panel: Fraction of LAEs in galaxies as a function of rest-frame UV (A = 1500 A) magnitude, predicted by the KTN09 model at z = 3. The 
dashed lines are the simple number fraction, while the solid lines are SFR-weighted fraction, i.e., the amount of star formation occurring in LAEs compared 
with the total star formation in galaxies. Thick and thin lines correspond to two different threshold values of rest-frame Lya EW of LAEs, 40 A and 120 A, 
respectively. Lower Panel: The rest-frame UV luminosity function of galaxies. The dashed line is the normal LF in number, but the solid and dotted lines are for 
the SFR-weighted and UV-luminosity- weighted LFs, respectively. The left and right ordinates are for the normal and SFR-weighted LF, respectively. 
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FIG. 4. — LAE fraction in GRB host galaxies as a function of rest-frame UV (A = 1500 A) magnitude of galaxies predicted by the KTN09 model at z = 3, for 
various values of critical metallicity Z cr j t (GRB occurs only when Z < Z cr i t ). The left and right panels are for the different threshold values of EWL, a for the 
LAE definition. Solid line is the case of no metallicity dependence, which is identical to the solid curves of the upper panel of Fig. [3] The 68 % C.L. confidence 
region from the observed data in Table[T]is shown by grey scale. 




FIG. 5. — LAE fraction in GRB host galaxies brighter than Muv -51og/i = — 21 (left) and —19 (right) in the KTN09 model at z = 3 is shown as a function of the 
model parameters (Z CI n or uz) of the metallicity dependence of GRB rate. The solid curve assumes that GRBs occur only when Z < Z cr i t , while the dashed curve 
assumes a Gaussian-type metallicity dependence of GRB rate with the dispersion parameter az- In the left panel, the 68% confidence region by the observed data 
in Table[T]is shown by grey scale. In the right panel, two models with different values of the threshold EWt ya are shown by thick (40 A) and thin (120 A) lines. 
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FIG. 6. — The same as Fig. [4] but here GRB afterglows are assumed to be detected only when Ay of their host galaxies is smaller than a critical value. The 
GRB rate is assumed to simply trace SFR without metallicity dependence. 



